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by  Brian  C.  McFall,  S.  Jarrell  Smith,  Cheryl  E.  Pollock, 
James  Rosati  III,  and  Katherine  E.  Brutsche 


PURPOSE:  This  Coastal  and  Hydraulics  Engineering  Teehnieal  Note  (CHETN)  deseribes  two 
methods  to  quickly  assess  sediment  mobility  relative  to  water  depth  as  a  funetion  of  loeal  waves 
and  eurrents.  Estimating  dredged  sediment  mobility  relative  to  placement  depth  and  grain  size  is 
valuable  for  preliminary  or  reeonnaissance  engineering  studies  to  compare  multiple  plaeement 
sites.  Presently,  a  Matlab  seript  is  available  to  perform  the  ealculations  using  both  methods.  Euture 
development  will  inelude  an  interactive  web  tool  with  automated  data  uploads  for  engineers  to 
utilize. 

BACKGROUND:  Nearshore  berms  are  eommonly  used  for  placement  of  dredged  sediment  that 
may  eontain  more  fine  silts  and  elays  than  are  allowed  for  placement  direetly  on  the  beach.  The 
United  States  Army  Corps  of  Engineers  (USACE)  produeed  several  Dredging  Researeh 
Technical  Notes  in  the  Dredging  Researeh  Program  (DRP)  Technical  Area  5,  Management  of 
Dredging  Projeets,  for  nearshore  berm  design  guidance  from  1990  to  1993.  The  majority  of  these 
teehnieal  notes  used  various  numerieal  models  for  design  guidanee.  A  brief  summary  of  the 
applieable  design  guidanee  obtained  from  the  teehnieal  notes  is  shown  in  Table  1. 


Table  1.  Summary  of  nearshore  berm  design  guidance  from  the  Dredging  Research 
Program  (DRP)  technical  notes. 

Technical  Note 

Key  Design  Guidance 

DRP-5-01 
(McLellan  1990) 

Recommended  a  linear  berm  rather  than  a  conical  berm  to  prevent  wave 
focusing. 

DRP-5-02 

(McLellan  et  al.  1990) 

Discussed  methods  to  determine  the  seaward  limit  of  the  littoral  zone  for 
active  berm  placement. 

DRP-5-06 

(Burke  and  Allison  1992) 

Recommended  minimum  longshore  length  of  610  m  (2,000  ft)  with  end 
slopes  of  1V:125H  to  prevent  wave  focusing. 

DRP-5-08 

(Pollock  and  Allison  1993) 

Recommended  minimum  berm  crest  width,  be,  61  m  (200  ft)  to  reduce  wave 
height. 

Although  the  praetiee  of  plaeing  dredged  sediment  in  the  nearshore  region  is  commonplace,  there 
is  not  a  simple  web  tool  applieation  available  that  ean  rapidly  produee  a  preliminary  assessment 
of  the  sediment  mobility  at  prospeetive  plaeement  sites  and  aid  with  the  evaluation  of  what 
portion  of  the  placed  sediment  is  likely  to  be  mobilized.  Estimating  sediment  mobility  is  the  first 
step  for  engineers  and  planners  to  estimate  the  potential  volume  of  sediment  that  a  placement 
operation  might  yield  to  beneficially  nourish  a  wetland  or  nearshore  region,  thereby  allowing 
stakeholders  to  gain  support  for  a  project. 
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METHODS  FOR  ESTIMATING  SEDIMENT  MOBILITY:  Two  methods  are  applied  to 
estimate  the  sediment  mobility  for  the  selected  location,  water  depth,  and  sediment  profile  of  the 
sediment  to  be  placed.  Method  1  analyzes  the  bed  shear  stress  from  local  wave  and  current 
conditions  and  compares  it  with  the  critical  thresholds  for  various  median  grain  size  diameters. 
The  bed  shear  stress  is  a  function  of  the  wave  orbital  velocity  and  is  calculated  using  linear  wave 
theory.  Method  2  analyzes  the  near-bed  velocity  and  compares  the  critical  near-bed  velocity  to 
locally  generated  velocity  data.  Method  2  uses  the  nonlinear  stream  function  wave  theory  to 
calculate  the  near-bed  wave  orbital  velocity,  which  generally  produces  larger  velocities  than 
linear  wave  theory.  Both  methods  are  applied  to  10  years  of  offshore  wave  hindcasts  from  the 
closest  Wave  Information  Studies  (WIS)  station  to  the  study  site  that  are  transformed  to  a  single 
nearshore  depth.  Additionally,  both  methods  are  applied  to  the  10-year  offshore  significant  wave 
height  and  period  data  that  is  then  transformed  to  a  range  of  nearshore  depths.  Both  Methods  1 
and  2  are  explained  within  this  technical  note  and  an  example  scenario  at  Duck,  NC,  is  provided. 
A  Matlab  script  is  available  to  automate  both  methods. 

Method  1  —  Bed  Shear  Stress.  This  method  imports  10  years  (1  January  1990-1  January 
2000)  of  wave  hindcast  data  from  the  closest  WIS  station  to  the  proposed  area  of  evaluation 
(http://wis. usace. army. mil/)  and  saves  it  as  a  .mat  file.  The  user  must  input  the  offshore  water 
depth  of  the  WIS  station,  the  depth  and  shoreline  orientation  of  the  nearshore  placement  location 
of  interest,  median  sediment  diameter  of  the  potential  nearshore  berm  sediment,  assumed  current 
velocity  1  meter  (m)  above  the  bed,  water  temperature,  and  salinity.  The  temperature  and  salinity 
data  are  used  in  the  density  and  viscosity  calculations. 

For  this  technical  note  an  example  study  site  is  selected  and  the  sediment  mobility  indexes  are 
calculated  using  the  local  shoreline  orientation  and  WIS  Station  63218  near  Duck,  NC.  The 
assumed  current  I  m  above  the  bed  is  O.I  meter  per  second  (m/s).  The  coordinate  system  of  the 
shoreline  angle  is  shown  in  Figure  1 . 


Figure  1 .  The  Atlantic  shoreline  orientation  at  Duck,  NC,  is  approximately  1 60°. 
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The  measured  wave  angle  and  height  from  the  WIS  station  are  transformed  from  offshore  to  the 
nearshore  region  using  eonservation  of  energy  flux  and  Snell’s  Law,  whieh  are  given  as 


and 


H,=H^ 


Cg^  |COS0„ 


V  COS0, 


■'g2 


(1) 


sin0^  _sin0„ 

respeetively,  where  H2  and  Ho  are  the  nearshore  and  offshore  wave  heights,  Cg  is  the  group 
veloeity,  0  is  the  angle  from  shore  normal,  C  is  the  wave  eelerity,  and  the  subseripts  “2”  and  “o” 
indieate  ealculations  in  the  nearshore  and  offshore,  respeetively.  The  polar  eoordinate  system 
used  for  9  is  shown  in  Figure  2. 


Figure  2.  The  polar  coordinate  system  used  in  the  sediment  transportation 
calculations. 


The  eritieal  shear  stress  is  estimated  from  Shields  diagram  following  a  proeedure  given  by 
Soulsby  (1997)  and  Soulsby  and  Whitehouse  (1997)  as 


1/3 


(3) 


0  =  —  +  0.055fl-exp(-0.020A 

1  +  1.2  A  ^  ^ 


(4) 


and 
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T^cr=^cr9{ps-p)ds0  (5) 

where  £),  is  the  dimensionless  grain  size,  dso  is  the  median  grain  size,  g  is  the  acceleration  due  to 
gravity,  is  the  sediment  density,  p  is  the  water  density,  v  is  the  kinematic  viscosity  of  water, 
is  the  Shields  parameter,  and  is  the  critical  shear  stress.  The  critical  shear  stress  is  the 

threshold  stress  for  which  the  sediment  can  be  expected  to  be  dislodged  from  the  seabed  for  all 
greater  shear  stresses. 

The  bottom  skin  shear  stress  is  calculated  using  a  method  described  by  Soulsby  (1997)  and 
Myrhaug  (1989)  for  currents  and  waves.  Form  shear  stress  is  not  included  in  the  calculations. 
The  current-induced  shear  stress,  ,  is  calculated  as 


Uk 


(6) 


where  U  is  the  assumed  mean  current  velocity,  k  is  the  von  Karman’s  constant  (k  =  0.4),  Zq  is 
the  bed  roughness  length  ( z^  =  / 12  for  flat  sand),  and  z  is  the  height  of  the  current  velocity 

assumed  by  the  user  (z  =  1  m).  The  wave-induced  shear  stress,  is  given  as 


T 


w 


(V) 


where  is  the  wave  friction  factor,  and  Uw  is  the  wave  orbital  velocity  determined  by 
Soulsby’ s  (1997)  orbital  bottom  velocity  figure  which  integrates  the  linear  wave  theory  orbital 
velocity  across  the  Joint  North  Sea  Wave  Observation  Project  (JONSWAP)  wave  spectra.  Linear 
wave  theory  is  appropriate  when  the  wave  steepness  (height/wavelength)  is  very  small  and  the 
magnitude  of  the  orbital  velocity  is  the  same  beneath  the  crest  and  the  trough.  The  friction  factor, 
f^,  is  calculated  using  the  Reynold’s  number,  Rw  ,  relative  roughness,  r,  and  the  semi-orbital 
excursion.  A,  which  are  given  as 


and 


R 


U...A 


U  T 

W 
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(8) 
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where  T  is  the  wave  period,  and  is  the  Nikuradse  equivalent  sand  grain  roughness 
{k^  =  30 Zq).  The  single  wave  period,  T,  is  applied  to  Equation  10  for  monoehromatic  waves, 
but  the  peak  wave  period,  Tp,  is  appropriate  for  spectral  wave  conditions  from  WIS  stations.  The 
friction  factor,  ,  is  calculated  as  the  greater  of  and  ,  which  are  given  as 

=0.237 r  (rough  turbulent  flow)  (11) 


and 


where 


=BRJ^  (laminar  and  smooth  turbulent  flow) 


(12) 


B  =  2,  N  =  0.5  for  <  5x10^  (laminar  flow)  (13) 


and 


B  =  0.0521,  N  =  0.187  for  R^  >  5x10^  (smooth  turbulent  flow)  (14) 


The  maximum  shear  stress,  ,  from  the  waves  and  currents  shear  stresses  is  calculated  as 


(  > 

3.2' 

1  +  1.2 

\  C  W  J 

(15) 


and 


T _ = 


T„,+T^cos(^f +(T^sin(^)^ 


1/2 


(16) 


where  is  the  mean  bed  shear  stress  and  c))  is  the  angle  between  the  wave  and  current 
directions. 

Figure  3  shows  the  histogram  of  maximum  bed  shear  stresses  using  transformed  wave  hindcasts 
from  WIS  Station  63218  to  a  nearshore  depth  of  8  m  near  Duck,  NC.  The  program  analyzed 
approximately  88,000  waves  and  required  112  seconds  (s)  to  run  after  the  WIS  data  were 
imported  into  a  .mat  file.  The  legend  shows  the  median  grain  size,  critical  shear  stress,  frequency 
of  mobility  during  the  10-year  period, /k,  and  mean  mobility  score  which  is  given  as 


M 


(17) 


The  frequency  of  mobility  applies  to  the  specified  median  grain  size  diameter.  Sediments  with 
smaller  grain  size  diameter  will  be  mobilized  more  frequently.  The  mean  mobility  score  allows 
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the  user  another  way  to  analyze  the  predietive  sediment  mobility  by  quickly  assessing  how  much 
the  maximum  bottom  stress  exceeds  the  critical  stress  on  average,  and  can  be  particularly  useful 
when  comparing  sites  with  similar  frequencies  of  mobility.  The  mean  mobility  score  can  be 
negative  in  sites  where  the  average  maximum  bottom  stress  is  less  than  the  critical  bottom  stress. 
Sites  with  negative  mobility  scores  will  experience  little  mobility. 

In  the  Duck,  NC,  example  shown  in  Figure  3,  the  critical  shear  stress  at  a  depth  of  8  m  ranged 
from  0.14  to  0.26  Pa,  resulting  in  71%  to  84%  of  the  waves  moving  the  sediment,  and  a  mean 
mobility  score  ranging  from  1.7  to  2.0.  The  finer  sediment  has  an  increased  frequency  of 
mobility  and  increased  mean  mobility  score,  as  expected.  For  this  example,  if  the  user  were 
placing  a  nearshore  berm  with  a  median  grain  size  of  0.2  millimeters  (mm),  it  could  be  estimated 
that  approximately  73%  of  the  waves  would  be  effective  in  transporting  sediment.  However,  a 
berm  with  median  grain  size  greater  than  0.5  mm  would  be  less  likely  to  be  transported. 


0  0.5  1  1.5  2  2.5  3 

Bed  Stress  [Pa] 

Figure  3.  Histogram  of  the  calculated  maximum  bed  shear  stress  at  a 
depth  of  8  m  with  the  critical  shear  stress  for  the  respective 
median  grain  sizes  noted  with  the  vertical  dashed  lines.  N  is 
the  number  of  waves  during  the  10-year  period  in  each  shear 
stress  bin. 

In  addition  to  transforming  tens  of  thousands  of  waves  to  a  single  nearshore  point,  the  offshore 
significant  wave  height,  period,  and  direction  are  calculated,  and  the  single  offshore  wave  is 
transformed  to  a  range  of  nearshore  depths  to  estimate  sediment  mobility.  The  range  of  nearshore 
depths  can  be  determined  by  the  user  based  on  the  depths  in  the  area  of  the  potential  berm  location. 

The  Matlab  code  imports  10  years  of  wave  measurements  from  the  closest  WIS  station,  calculates 
the  significant  wave  height  and  period,  and  applies  Equations  1-16  to  calculate  the  maximum  bed 
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shear  stress  for  the  range  of  nearshore  depths.  Figure  4  shows  the  maximum  bed  shear  stress  indueed 
by  the  offshore  signifieant  wave  height  Ho  =  1.82  m,  peak  offshore  wave  period  Tp  =  8.89  s,  and 
shoreline  angle  of  160°.  The  maximum  bed  shear  stress  indueed  by  this  wave  height  and  period 
exceeded  the  critical  thresholds  for  mobility  at  all  nearshore  depths  from  3  to  12  m. 


H  =1.82m,T  =  8.98  s 

o  p 


Bed  Stress  [Pa] 


Figure  4.  The  total  maximum  bed  shear  stress  for  nearshore  depths  from 
3  to  12  m  for  the  significant  offshore  wave  height  and  period 
from  WIS  Station  63218. 

Method  2  —  Near-Bottom  Velocity.  This  method  imports  10  years  (1  January  1990-1 
January  2000)  of  wave  data  from  a  WIS  station,  but  determines  the  sediment  mobility  using  the 
near-bottom  critical  velocity  instead  of  the  bed  shear  stress.  The  near-bottom  velocity  is 
calculated  using  nonlinear,  stream  function  wave  theory.  By  also  calculating  the  sediment 
mobility  using  this  method  as  well  as  Method  1,  a  range  of  predicted  mobility  is  provided. 
Stream  function  wave  theory  is  appropriate  for  water  depths  in  the  range  0.006  gP  to  0.016  gP 
(Soulsby,  1997),  or  for  the  Duck,  NC,  example  in  the  range  4.7  to  12.6  m. 

The  critical  near-bottom  velocity,  Ucr  as  given  by  Ahrens  and  Hands  (1998)  is  based  on  research 
by  Hallermeier  (1980)  and  Komar  and  Miller  (1974),  is  given  as 

Ucr  =^J^9Yd^o  d^o  <  2.0mm  (18) 


and 


u 


cr 


4/7 


for  dgQ  >  2.0mm 


(19) 
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where  y  is  defined  as  y  =  (p^  —  p)/p ,  where  p^.  is  the  sediment  density  and  p  is  the  water 

density.  Ahrens  and  Hands  (1998)  used  Dean’s  (1974)  stream  funetion  wave  theory  table 
(SFWT)  to  derive  the  following  equations  for  the  near-bottom  wave  indueed  veloeity  based  on 
stream  function  wave  theory  for  the  wave  crest,  and  trough, 


and 


H 

f/i] 

^maxcrest 

T) 

-0.579 


exp 


0.289-0.491 


H 

li 


2.97 


u 


max  trough 


H 


exp 


1.966-6.70 


-1.73 


H 


+  5.58 


H 


(20) 


(21) 


where  h  is  the  water  depth,  H  is  the  wave  height  in  the  placement  site,  and  is  the  offshore 
wave  length  given  by  L^  =  {gT^^/27i.  The  maximum  near-bottom  velocity  was  taken  as 


=  max  u 


max  trough 


)■ 


Figure  5  shows  a  histogram  of  the  near-bottom  velocity  as  a  function  of  the  number  of  waves  in 
each  shear  stress  bin.  The  critical  velocities  per  median  grain  size  are  shown  with  dashed  lines. 
The  program  required  6  seconds  to  analyze  88,000  waves  after  the  WIS  data  were  imported  into 
a  .mat  file.  The  frequency  of  mobility  ranged  from  84.5%  to  98.7%  depending  on  the  median 
grain  size.  The  mean  mobility  score  using  the  near-bottom  velocity  is  calculated  as 


= 


Wmax-Wc: 


u„ 


(22) 


Figure  5  is  a  great  example  of  the  importance  of  the  mean  mobility  score.  For  median  grain  sizes 
of  0.1  and  0.2  mm,  the  frequency  of  mobility  is  within  3%,  but  the  mean  mobility  scores  are 
significantly  different  (2.6  vs.  4.1).  The  greater  mobility  score  indicates  that  the  smaller  grain 
size  will  move  farther  when  mobilized  and  most  likely  winnow  the  finer  sediment  from  the  site. 

The  near-bottom  maximum  and  critical  velocities  (Equations  18-21)  were  applied  to  multiple 
depth  contours  with  the  single  offshore  significant  wave  height  and  period  calculated  from  10 
years  of  offshore  WIS  wave  measurements.  Figure  6  shows  the  near-bottom  velocities  for  a 
range  of  nearshore  depths.  The  critical  velocities  for  several  median  grain  sizes  are  shown  with 
the  vertical  dashed  lines.  The  near-bottom  velocities  for  the  range  of  depths  shown  exceed  the 
critical  velocity  for  all  median  grain  sizes,  so  the  user  would  anticipate  that  berms  with  these 
median  diameters  would  be  mobilized. 


8 


ERDC/CHL  CHETN-IV-108 
April  2016 


-  -  d„  =  0,1  mm,  u  =0.11  m/s.  f,.  =  98.7  %,  <M 

50  cr  Mu  u 

Mu 


—  —  d„  =  0.2  mm,  u  =  0. 16  m/s,  f,.  =  96.0  %,  <M 

50  cr  Mu  u 

—  -d„  =  0.3mm,  u  =  0.19  m/s,  f,.  =92.3%,  <M 

50  cr  Mu  u 

—  —  d„  =  0.4  mm,  u  =  0.22  m/s.  f, .  =  88.3  %,  <M 

50  cr  Mu  u 

- d„  =  0.5  mm,  u  =  0.25  m/s,  =  84.5  %,  <M 

50  cr  Mu  u 


>  =  4.124 

>  =  2.623 
>=  1.958 
>=  1.562 
>=  1.292 


0  0.5  1  1.5  2  2.5  3 

Near-Bottom  Velocity  [nrVs] 


Figure  5.  Histogram  of  the  calculated  maximum  near-bottom  velocity  with 
the  critical  velocity  for  the  respective  median  grain  size  noted 
with  the  vertical  dashed  lines. 
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50  cr 
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18 


Near-Bottom  Velocity,  u^  [m/s] 


Figure  6.  The  total  maximum  near-bottom  velocities  for  nearshore 
depths  from  3  to  12  m  from  the  offshore  significant  wave 
height  and  period.  The  critical  velocity  for  each  median  grain 
size  is  shown  with  vertical  dashed  lines. 
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SUMMARY:  This  technical  note  discusses  two  methods  to  estimate  sediment  mobility.  Both 
methods  are  applied  to  10  years  of  offshore  wave  hindcasts  which  are  transformed  to  a  single 
nearshore  depth,  and  the  offshore  significant  wave  height  and  period  are  transformed  to  a  range  of 
nearshore  depths.  Figures  are  created  that  show  the  critical  shear  stress  and  critical  near-bed 
velocity  for  a  range  of  median  grain  sizes,  the  frequency  of  sediment  mobility  with  relation  to  grain 
sizes,  and  a  general  mobility  index.  Method  1  analyzes  the  bed  shear  stress  and  Method  2  analyzes 
the  near-bed  velocity.  The  bed  shear  stress  analyzed  in  Method  1  is  a  function  of  the  near-bed 
velocity.  The  near-bed  wave  orbital  velocity  for  Method  1  is  calculated  using  linear  wave  theory, 
which  is  appropriate  when  the  wave  steepness  (height/wave  length)  is  small.  Method  2  calculated 
the  near-bed  orbital  wave  velocity  using  stream  function  wave  theory,  which  is  appropriate  for 
water  depths  in  the  range  0.006  gT^  to  0.016  gT^,  or  for  the  Duck,  NC,  example  in  the  range  4.7  to 
12.6  m.  The  orbital  wave  velocity  is  larger  using  the  nonlinear  stream  function  wave  theory 
compared  to  linear  wave  theory,  hence  the  85%  to  99%  frequency  of  mobility  shown  in  Figure  5 
compared  to  the  71%  to  84%  frequency  of  mobility  shown  in  Figure  3.  Waves  become  more 
nonlinear  as  they  transform  into  the  nearshore,  making  Method  2  more  appropriate  and  Method  1 
more  conservative  for  the  example  herein.  By  using  both  methods,  a  range  of  mobility  is  shown. 

FUTURE  DEVELOPMENT:  The  two  methods  described  in  this  technical  note  will  be 
developed  into  a  web  application  to  provide  engineers  a  rapid  reconnaissance  tool  for  the 
preliminary  assessment  of  nearshore  berm  locations.  If  these  methods  are  applied  to  a  web  tool, 
the  user  would  select  the  WIS  station,  and  a  default  value  for  the  temperature,  salinity,  and 
current  1  m  above  the  bed  would  appear  for  the  user.  The  web  tool  would  include  direct  data 
uploads  for  the  selected  WIS  station.  The  depth  of  closure  calculation  could  be  included  in  this 
tool  using  the  method  and  program  described  by  Brutsche  et  al.  (in  prep.).  The  information 
provided  from  the  tool  could  potentially  be  used  to  create  a  map  with  colored  contours  to  show 
the  sediment  mobility  at  various  depths. 

ADDITIONAL  INFORMATION:  This  Coastal  and  Hydraulics  Engineering  Technical  Note 
(CHETN)  was  prepared  as  part  of  the  USAGE  Coastal  Inlets  Research  Program  (CIRP)  by 
Dr.  Brian  C.  McEall,  Dr.  S.  Jarrell  Smith,  Cheryl  E.  Pollock,  James  Rosati,  III,  and  Dr.  Katherine 
E.  Brutsche,  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC),  Coastal  and 
Hydraulics  Eaboratory  (CHE),  Vicksburg,  MS.  Questions  pertaining  to  this  CHETN  may  be 
directed  to  Brian  McEall  (Brian.  C.McFall(a),usace. army. mil)  or  to  the  USAGE  CIRP  Program 
Manager,  Julie  Dean  Rosati  ( Julie. D.Rosati(a), usace.armv.mil).  Additional  information  regarding 
CIRP  may  be  obtained  from  the  CIRP  web  site  http://cirp. usace. armv.mil/. 

This  ERDC/CHE  CHETN-IV-108  should  be  cited  as  follows: 

McEall,  B.  C.,  S.  J.  Smith,  C.  E.  Pollock,  J.  Rosati,  III,  and  K.  E.  Brutsche.  2016. 

Evaluating  sediment  mobility  for  siting  nearshore  berms.  ERDC/CHE  CHETN- 

IV-108.  Vicksburg,  MS:  U.S.  Army  Engineer  Research  and  Development  Center. 
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